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ABSTRACT 

We present the mass-X-ray observable scaling relationships fo r clusters of 



galaxies using the XMM-Newton cluster catalog of iSnowden et all Our results 
are roughly consistent with previous observational and theoretical work, with 
one major exception. We find 2-3 times the scatter around the best fit mass 
scaling relationships as expected from cluster simulations or seen in other ob- 
servational studies. We suggest that this is a consequence of using hydrostatic 
mass, as opposed to virial mass, and is due to the explicit dependence of the 
hydrostatic mass on the gradients of the temperature and gas density profiles. 
We find a larger range of slope in the cluster temperature profiles at r 500 than 
previous observational studies. Additionally, we find only a weak dependence of 
the gas mass fraction on cluster mass, consistent with a constant. Our average 
gas mass fraction results argue for a closer study of the systematic errors due 
to instrumental calibration and analysis method variations. We suggest that a 
more careful study of the differences between various observational results and 
with cluster simulations is needed to understand sources of bias and scatter in 
cosmological studies of galaxy clusters. 

Subject headings: cosmology: observations — galaxies: clusters: general — X- 
rays: galaxies: clusters 
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Introduction 



Studies of clusters of galaxies provide for a variety of cosmological tests (see IVoitl 120051 . 
for a recent review). The precision of these tests is limited by the accuracy and precision of 
the scaling relations used to transform X-ray observables, such as temperature or luminosity, 



into cluster mass measurements (e.g., iRapetti et al.ll2008l ; IVikhlinin et al.l 120091 ) . Improved 



modeling of the complex physics at work in clusters has led to better study of the sources of 
bias and scatter in the mass scaling relationships. These models suggest that the hydrostatic 
equilibrium assumption used to calculate total c luster masses from X-ray data can lead 



to underestimates on the order of 10-20% (e.g., iNagai et al.l l2007bl ). The discrepancy is 
explained by the presence of nonthermal pressure support in clusters. Theoretical studies 
have found that the expected scatter in mass sc aling relationships is dependent on the 



dynamical state of the clusters under study (e.g., iKravtsov et al.l 120061 ) . Relaxed systems 



should show lower scatter around the best-fit relation than disturbed clusters. 

This information can be used to tailor observational programs, particularly by limiting 
the bias and scatter in the mass scaling relationships. One suggestion is to use only relaxed 
clusters in cosmological studies. However for high r edshift samples , relax ed clusters make up 
only a minority of the cluster population (see e.g., IVikhlinin et al.l 120091). Altern atively, the 
right choice of mass proxy could provide low scatter data. IKravtsov et al.l (120061 ) proposed a 
new X-ray proxy for cluster mass, the parameter Yx which is the product of X-ray temper- 
ature and gas mass. The Yx parameter is related to the total thermal energy of the cluster 
gas and is an X-ray a nalog to the Sunyaev-Zel'dovich (SZ) flux. From cluster simulations, 
Kravtsov et al.l (120061 ) found that the scatter in the cluster mass-Yx scaling law is not only 
lower than for other commonly used mass proxies, but also shows little dependence on cluster 
dynamical state. 

Observational studies of the Yx parameter as a mass proxy have been limited in the num 



berof clusters used and the statistical quality of the data (lArnaud et al.ll2007l ; IVikhlinin et al. 



20091 ). The statistical error was on the order of the measured scatter, making it difficult 
to determine the intrinsic scatter in the relationship. Additionally, previous observational 
work has focused solely on a limited number of relaxed systems. A larger and more sen- 
sitive observational study of clusters will better test the Yx mass proxy and the results of 
cluster simulations. In this Letter, we use the data from a recent high signal-to-noise XMM- 
Newton survey of nearby clusters to test the usefulness of various mass proxies. We assume 
H = 70 km s" 1 Mpc" 1 , Q M = 0.3, and Q A = 0.7. 
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2. Data Analysis 



Our cluster sample was taken from lSnowden et al.l (120081 ). which presented the projected 
temperature, abundance, and surface brightness profiles for 70 clusters found by fitting 
XMM-Newton data (see that work for details of the spectral analysis). While not an unbiased 
sample, the large number of clusters does provide a wide range of cluster properties to 
study. I t is not limited to re laxed and/or hot (> 5 keV) clusters. The selection criteria 
used by ISnowden et al.l (120081 ) excluded highly asymmetric clusters and those with strong 
substructure (e.g., A115, A754). However, some merging clusters were included in the sample 
(e.g., A520, the Bullet Cluster). 

We used the projected profiles to deter mine the gas dens i ty an d three-dimensional tem- 
perature profiles following the procedure of IVikhlinin et al.l (120061 ). The same gas density 
and temperature models were used with some simplification. The lower spatial resolution of 
our data did not constrain the second /3-model component required to fit the Chandra data 
in the cluster center, therefore we did not include it. We also compared fits with and without 
steepening at large radii and excluded the steepening when it did not produce a significantly 
better fit of the data. For the temperature fits, we t ook into account projection effects and 
used the spectral temperature weighting formalism of lVikhlininl (120061 ). When cooling in the 
core was not obvious we fixed the cool component parameters to produce a flat temperature 
profile in the center. 

We calculated the total cluster mass distribution for each cluster assuming hydrostatic 
equilibrium and the radii (r 2 5oo and r 500 ) where the cluster density equaled 2500 and 500 
times the critical density at the cluster redshift. We do not include clusters where the 
calculated r 2 soo and r 500 values extend beyond the radial coverage of our data. We find that 
60/70 clusters have data extending to at least r 2 soo and 28/70 have data out to r 50 o. We 
determined the gas mass (M g ) and total cluster mass (M) enclosed by r 2 5oo and r5oo for each 
cluster. We also calculated the aver age spectral tem perature (T x ) within the 0.15-1 r 500 
radial range using the formulation of I Vikhlininl (120061 ). 



Uncertainty intervals were obtained from Monte Carlo simulations. We simulated sur- 
face brightness and projected temperature distribut ions by scatter i ng th e observed data 
according to the measurement uncertainties found in ISnowden et al.l (120081 ) . The simulated 
data were fit with gas density and temperature models and a full analysis performed to 
determine M, M g , and Tx- The uncertainties were obtained from their distribution in the 
simulated data. For values evaluated at r 50 o (^2500), the uncertainty includes the uncertainty 
on r 50 o (r 25 oo)- 
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3. Comparison of Scaling Relations 



We determined the best-fit scaling relationships at r 50 o using E(z) n M = C(X/X ) a for 
X = Tx, M g , and Yx- We fixed n — 1, and 2/5, as consistent with our cosmology, and 
X = 5 keV, 4 x 10 13 M and 3 x 10 14 M keV for the T x , M g and Y x fits, respectively. We 
also present the best-fit relationship between the gas mass fractio n, f g = M g /M, and M , as 



sec 



Vikhlinin et all 120091 ). 



characterized by the equation f g = / S)0 + alog 10 [M/10 M 

We use the BCES fitting routines^ which provide a linear regression alg orithm that al- 
lows for intrinsic scatter and nonuniform measurement errors in both variables (lAkritas fc Bershady 
19961 ). We find that the Y|X and orthogonal slope estimators provide the most significant 
(and consistent) results compared to other methods, including bisector. We therefore only 
present the results of the Y|X and orthogonal methods for each of our fits (see Tabled]). We 
also include the best-fit results when the relationship slope, a, is fixed at the expected value 
from self-similarity. 

For the M — X relationships we esti mate the intrins i c scat ter using a generalized form 
of the estimated scatter, 5M/M, used by I Vikhlinin et al.l (120061 ): 



5M 



N-2 



[Mj - C^/XqY} 2 - AM? 
M? 



where AMj are the measurement errors. Similarly for the f g — M relationship, we calculate 
the scatter Sf g / f g by: 

2 



6A 



i 



N-2 



\- [fg,i - (fg,0 + alog lo [M t /lO 15 M ])] 2 - Afl 



where Af g ^ ar e the measurement errors. To compa re our scatter with logarithmic scatter 
estimates (e.g.. |Jeltema et al.ll2008l ; iPratt et al.ll2009l ). multiply the logarithmic estimates by 
In 10 = 2.30. Given the high statistical quality of our data, the scatter is dominated by the 
intrinsic scatter. 



3.1. Mass— Temperature 



Table [T] and Figure [T] show the best-fit parameters for the M — Tx scaling relation. 
Our best fit (l og in C = 14.64, a = 1.67) is consistent w ith the scaling relations found in 
other datasets fkrnaud et ali 120071 : 1 Vikhlinin etaJbood . log 10 C = 14.580 and 14.635, and 



Routines available at 



http://www.astro.wisc.edu/~mab/archivc7stats/stats.html 
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a = 1.71 and 1.53, respectively). lArnaud et al.l (120071 ) use a different definition of Tx which 
integrates the observed temperature profile ov er 0.15-0.75 r^nn- Since t he average cluster 
temperature profile falls at large radii (see e.g.. iLeccardi fc Molendill2008l . A. M. Juett et al. 
2009, in prep), this definition will produce higher values of T x and subsequently lower values 
of the M -T x normalization. We find that (T x (0.15 - 0.75r 5 oo)/T x (0.15 - lr 500 )) = 1.06, 
which translates into a reduction of lo g in C o f 0.04-0.05, enough to explain the discrepancy 
between our results and lArnaud et al.l (120071 ) . 

When compare d to theoretical calculations ( Kravtsov et al. 2006 ; Nagai et al. 2007al lbl; 



Jeltema et al.l 120081 ). our results are consistent when the hydrostatic mass is considered, 
but our normalization is lower by « 20% when compared to scaling results that use the 
true cluster mass (14.70-14.75). This is a well known result and is likely due to non-thermal 



pressu re support that is not accounted for in the hydrostatic mass estimate (e.g. jNagai et al. 



2007al ). The biggest difference between our results and previous studies is the difference 



in scatter around the best-fit scaling relation. We find a scatter of 43% which, due to 
our systematics limited dataset, can be attributed to the intrinsic scatter of our sample. 
Other studies, of both observed and theoretical cluster samples , have found significantly 
lower values for th e i ntrinsic scatter, pa 1 — 25% (jVikhlinin et al.ll2006l ; lArnaud et al.l 120071 ; 
Nagai et al.ll2007al lbl: [jeltema et al.lboosf ). 



3.2. Mass- Gas Mass 



The gas mass-total mass scaling relationship has been previously characterized in two 
ways: a power law scaling between M and M g , and more recently by a linear scaling of f g 
and the logarithm of M. 

Our M — M g results (log_m C = 14.503±0.018), are in agreement with the normalization 



found by lArnaud et al.l (120071 . 14.542±0.015) considering the small formal errors. But our 
best-fit slope is steeper (0.93±0.05 vs 0.80±0.04). Again we find that theoretical models 
that use the true cluster mass have higher predict ed normalizations but our results are con - 
sistent when hydrostatic mass estimates are used (IKravtsov et al.ll2006l ; iNagai et al.ll2007al ). 
The slope estimate is consistent with theoretical results for both true and hydrostatic mass 
estimates. Our 15% scatter is close to the ~10% scatter found in both observational and 
theoretical work. Interestingly, our res ults match t h e comb ination of normalization, slope, 
and scatter found in the simulations of INagai et al.l (j2007al ) when hydrostatic mass is used, 
however we differ from their Chandra observational results. Our slope is larger (0.93 vs 0.81) 
while our normalization is lower (14.503 vs 14.59). This discrepancy may be d ue to differences 



in the calibration of the instruments that has been previously noted (e.g., ISnowden et al. 
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2008|). 



We find that the normal i zation of the f g —M relationship, 0.134±0.005, is consistent with 
the value of IVikhlinin et al.l (120091 . 0.130±0.007). Our data are con s istent with a constant 
slope over the range of mass considered, w hile the IVikhlinin et al.l (120091 ) result prefers a 
reduction in f g for lower mass clusters. The IVikhlinin et al.l (120091) res ult is consistent with 
gas mass fraction results from groups of galaxies (jSun et al.l 120091 ) . lArnaud et al.l (120071 ) 
suggested that the gas mass fraction may be constant above 2 — 3 x 10 14 M Q , and then 
dropping at lower masses. This would explain both our result and the lower gas mass 
fractions seen in groups of galaxies. 



The differe nce in mean gas mass fraction between the IVikhlinin et al.l (120061 ) work, 
a subset of the IVikhlinin et al.l (120091 ) sample, and ours is 10-20%. At r ^ nn, w e find a 
mean f g = 0.1323 ± 0.0019, compared to 0.110±0.002 for the IVikhlinin et all tood ) sample 
(adjusted to account for differences in cosmology). If we restrict our analysis to clusters with 
kT > 5 keV, the diff erence is «10%, 0.138±0.003 for our work compared to 0.123±0.003 
Jvikhlinin et alil2006h . 



Our results are also close to those found by lAUen et al.l (120081 ). At r 2 5oo, lAUen et al. 



(120081 ) found a mean cluster gas mass fraction of 0.113±0.003 for clusters with kT > 5 keV 
and low redshifts (z < 0.15). For their full sample (k T > 5 keV but all redshifts) they find 
a mean f g = 0.1104 ± 0.0016. IVikhlinin et al.l (120061 ) noted that their mean f g at r 250 o was 
significantly less (0.091±0.002, a ~25% difference) than an earlier (but consistent) Allen et 
al. sampl e. We find f n^noa — 0.1057 ± 0.0005 for clusters with kT > 5 keV, a 4% dif f erence 
with the lAUen et al.l (120081 ) results and a ~15% difference with the IVikhlinin et al.l (120061 ) 
work. 



3.3. Mass-F; 



x 



Kravtsov et al.l (120061 ) suggested that a lower scatter (< 10%) proxy for cluster mass 
is the parameter Yx = TxM g . Table [1] gives our best-fits for the M — Yx sc aling relation. 
Our b est-fit normalization (14.657±0.018) is consistent with the r esults of lArnaud et al. 
J2007L 14.653±0.015 when renormalized) and IVikhlinin et all (I2009I . 14.684±0.015). When 
comparing with theoretical results, we find consistency wit h hydrostatic mass results (14.645) 
but not true cluster mass (14.712; e.g. iNagai et al.ll2007al ). The slope of the best-fit M — Y x 
compares well with both observational and theoretical studies. 

The largest difference comes in the measured scatter in the M — Yx relation. We find 
a scatter of ~ 22%, much larger than the < 10% level expected from simulations. While 
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other observational studies have not found as large a scatter (e.g., lArnaud et al.l 120071 ). we 
note that ours is the first study to be systematics limited and includes twice the number of 
objects. Thus we are able to measure the scatter without a significant contribution from the 
statistical error. In theoretical wo rk, scatter increases to 8-20% when hy drostatic masses 
are used compared to true masses (jNagai et al.ll2007al ; IJeltema et al.ll2008l ). 



4. Correlation of Deviation with Cluster Properties 



Scatter in the mass-X-ray observable relationships is an important contributo r to the 
total error budget in cosmological studies of clusters (see e.g. IVikhlinin et al.l 120091 ). Given 
our large sample of clusters, we can study what factors are most important in producing 
scatter in these relationships which can then be used to refine cosmological studies to reduce 
the scatter. 

One possible cause of the observed scatter is the dynamical state of the cluster. Relaxed 
clusters are expected to better follow the assumptions of hydrostatic equilibrium. Disturbed 
clusters may have additional pressure and energy input s due to the mergin g events and the 



additional complication of asymmetric geometries (e.g.. lNagai et al.ll2007bl ). 



We looked for correlation between the deviation of the calculated hydrostatic mass, 
Mi, from the expected mass given the best-fit scaling relationship, M(Xi), with measures 
of the ellipticity and asymmetry in the cluster images (see Figure [2]). We identify the 
deviation as 8M(X) / M = \Mj — M(Xj) ] /Mj. The ellipticity and asymmetry were calculated 
following the work of lHashimoto et aJj (120071 ). Relaxed clusters should have low ellipticity 
and asy mmetry values, w h ile dis turbed systems will have higher ellipticity and/or asymmetry 
values. lHashimoto et al.l (120071 ) showed that ellipticity is correlated wi th the P2/P0 power 



ratio and that asym metry is related to the P3/P0 power ratio (see e.g., iBuote h Tsail 11995 



PC 

1 



Jeltema et al.ll2008l . for a discussion of power ratios). We find no correlation between the 
amount of deviation and either ellipticity or asymmetry for any of our scaling laws. 



We looked for other cluster properties that might influence the deviation. IVikhlinin et al. 



( 120061 ) noted that implicit in the calculation of hydrostatic mass, there is a dependence 
of the normalization of the M — Tx relation on the sum of the temperature gradient, 
fit = (—1/3) d log T/d log r, and gas density gradient, (3 e ff = (—1/3) c/logp/c?logr (see their 
Appendix A). If the spread in values is large enough, this dependence should cause a pre- 
dictable deviation around the best-fit relationship. We find a strong correlation between 
5M(X)/M and /3 e jf + fit for all scaling relationships (Figure EJ). The Spearman rank cor- 
relation coefficients were 0.62, 0.53, and 0.61 for X = T x , M g , and Y x , corresponding to 
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probabilities of 0.9987, 0.9943, and 0.9986, respectively." 



The I Vikhlinin et al.l (120061 ) cluster sample has a narrow distribution of /3 e ff + j3 t values. 
Our sample however, shows a large variation in the temperature gradient. For (3 t , we find 
a mean of 0.53 and a standard deviation of 0.50. The range of temperature gradient values 
i s a reflection of the variation of temperature profiles at large radii in our cluster sample 
( jSnowden et al.l 120081 ) . The gas density gradient has a narrower distribution with a mean 
Peff = 0.63 and a standard deviation of 0.10. Our gas density gradient results are in 
good agreement wi th those found in the REXCESS study (f3 ne (0.3-0.8 r 500 ) = 0.60±0.10 
Croston et aDl2008h . 



5. Discussion 



Our XMM-Newton survey of galaxy clusters (jSnowden et al.ll2008l ). provides a large 
sample and high signal-to-noise data to study the scaling relationships between cluster mass 
and cluster X-ray temperature, gas mass, gas mass fraction, and the mass proxy Yx- Our 
fits are in good agreement with other observational and theoretical work, with one major 
caveat. We find a significantly larger scatter around the best- fit relationships than has been 
previously seen. 

The scatter around the best-fit scaling relationships is 2-3 times that found in most 
other observational and theoretical work. Interestingly, a recent study of the scaling rela- 
tionship between the Compton y-parameter from SZ studies and cluster masses obtained 
from gravitational lensing also showe d large scatter in th e Mql ~ Tx and Mql ~ Y rela- 
tionships (41% and 32%, respectively; IMarrone et al.ll2009T ) , how ever we note a difference in 
methodology between their work and ours (IMarrone et al.l (120091 ) used a fixed clustercentric 
radius of 350kpc, while we work in mass scaled units). While other w ork has suggested tha t 



disturbed clusters will show more scatter than relaxed clusters (e.g. iKravtsov et al.ll2006l ). 
those authors used a visual classification of their simulated clusters, rather than a quantita- 
tive determination, making comparison difficult. Within our sample, we find no correlation 
between measures of the cluster dynamical state (ellipticity and asymmetry) and deviation 
around the best-fit scaling relationships. 

Only one cluster property, the combination /3 e tf + (3 t , showed a strong correlation with 
deviation. Since hydrostatic mass estimates, like those used here, depend on this explicitly, 



the result should be expected (e.g. IVikhlinin et al.ll2006l ). Our data are the first to show this 



dependence due to the large sample size and range of cluster properties included. There is 
no reason to suspect that true cluster masses will show such a dependence given the bias 



- 9- 



(and scatter) in hydrostati c mass determinatio ns found in simulations where the true cluster 
masses are known (see e.g. iNagai et al.ll2007bl ). 



One question we must ask is how reliable are our determinations of the temperature 
profiles, whose wide range dominates the measurement of P ff + t - To check the reliability 
of the background modeling and spectral fitting procedure, ISnowden et al.l (120081 ) compared 
the profile of A1795 with published Chandra, XMM-Newton, and Suzaku temperature profiles 
and found no significant difference. They also found reproducibility for three clusters with 
multiple observations. An initial comparison of our average cluster temperature profile is 
in good agreement with previous results both in overall shape and expec ted scatter (see A. 
M. Juett et al. 2009, in prep). The results of iLeccardi &: Molendil (120081 ) also suggest that 
cluster temperature profiles, while generally showing a falling profile at large radii, do show 
a range of profiles. Followup observations are needed of the most unusual systems to confirm 
our results. 

Assuming our temperature profile range is indicative of the cluster population, we then 
need to ask how does this result affect cosmological studies using clusters. First, a larger 
scatter should b e taken into considerat ion when discussing systematic errors in cosmological 
studies (see e.g. I Vikhlinin et al.l 120091 ) . However, given other error sources, it is not clear 



that the mass scaling relationship scatter would be the dominant error contributor. 

It may be possible to correct for the scatter from theoretical st udies of clust e r prop erties 
but it is unclear if present models are consistent with our results. INagai et al.l (j2007af ) find 
little variation in the temperature profile at r 50 o in their simulations, although these are 
limited to their relaxed subsample. If a more thorough study of the simulations is not able 
to reproduce the observed cluster variation, that may point to some missing physics needed 
to better describe the conditions within clusters. 

Another result we would like to highlight is the comparison of our averag e gas mass frac 



tion with previous resu lts. The differences between our work and others (e.g.. I Vikhlinin et al. 



20061 ; lAllen et al.ll2008l ) range from 5-20%. This is comparable with expected systematic dif- 
ferences between the instruments and analysis methods, but is significantly larger than the 
statistical errors typically quoted. In our opinion, a study of the expected systematics due to 
(1) instrumental calibration differences, and (2) data analysis methods must be performed. 
These issues are beyond the scope of this work, but will be addressed in our future paper 
(A. M. Juett et al. 2009, in prep). 
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Fig. 1. — Upper Left: Plot of X-ray spectral temperature, Tx, and total cluster mass, M. 
Overplotted are the best-fit power-law relations using the BCES orthogonal slope estimator 
(red), the BCES Y|X slope estimator (blue), and a fixed slope of a = 1.5 (green). Upper 
Right: Plot of Yx and M with best-fit power-law relations overplotted. Color coding is the 
same as for the M — Tx plot with a = 0.6. Lower Left: Plot of cluster gas mass, M g , and 
M and best-fit power-law relations overplotted. Color coding is the same as for the M — Tx 
plot with a = 1.0. Lower Right: Plot of gas mass fraction, f g , and M. Overplotted is the 
best-fit relation. 
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Table 1. Best-Fit Parameters of Powerlaw Fits 



Relation a 


Fit Method b 


!Ogl C/fgfl 


a 


Scatter 


M-T 


Y|X 


14.64±0.03 


1.67±0.16 


0.420 




Orth 


14.63±0.03 


1.89±0.20 


0.430 




Fix 


14.63±0.03 


1.5 


0.412 


M-Y x 


Y|X 


14.657±0.018 


0.61±0.04 


0.218 




Orth 


14.657±0.018 


0.61±0.04 


0.218 




Fix 


14.646±0.018 


0.60 


0.212 


M-M 9 


Y|X 


14.503±0.018 


0.93±0.05 


0.144 




Orth 


14.504±0.017 


0.93±0.05 


0.145 




Fix 


14.484±0.013 


1.0 


0.149 




Y|X 


0.134±0.005 


0.011±0.013 


0.104 




Orth 


0.134±0.005 


0.011±0.013 


0.104 



a The relationships are of the form E(z) n M = C(X/X ) a for the M-T, 
M - Y x , and M - M g fits with n = 1, 2/5, and 0, and X = 5keV, 3 x 
10 14 M o keV and 4 x 10 13 M & for the T, Y x , and M g fits, respectively. For 
the / fl -M fit, the relationship is f g = f gfi + a log lo (M/lO 15 M ). 

b We use the BCES fitting package and present results from the Y|X and 
orthogonal fitting methods. In addition, we fit the data with a fixed slope, 
a, given by the expected self-similar relationships (a = 1.5, 1.0, and 0.6 for 
the T, Y x , and M g fits, respectively). 



>- 
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Fig. 2. — Plot of deviation (5M(Y X )/M) around the best-fit M — Y x scaling relationship 
versus measures of the cluster dynamical state, ellipticity (left panel) and asymmetry (right 
panel). No correlation is found. 
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Fig. 3. — Plot of deviation around the best-fit mass scaling relations, 8M(X)/M, versus 
fleff + A- The results for the scaling relationships for M — Tx (top panel), M — M g (middle 
panel) and M — Y x (bottom panel) are given. All show a correlation between deviation and 

Peff + Pf 



